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Abstract — In this paper, the comparative simulation results 
of using a static synchronous compensator (STATCOM) and a 
static VAR compensator (SVC) to achieve damping 
enhancement of a three-machine system is presented. PID 
damping controllers of the proposed STATCOM and SVC are 
designed by using Modal Theory Method to contribute adequate 
damping characteristics to the dominant modes of the studied 
system. A frequency-domain approach using eigenvalue analysis 
and a time-domain scheme based on a nonlinear system model 
subject to a three-phase short-circuit fault at the load connected 
bus are systematically performed to examine the effectiveness of 
the proposed control schemes. It can be concluded from the 
comparative simulated results that the designed PID damping 
controller joined with the proposed STATCOM is shown better 
than the proposed SVC. 

Index Terms — Static Synchronous Compensator 
(STATCOM), static VAR compensator, Multi-machine System, 
Stability. 


I. INTRODUCTION 

Nowadays, flexible AC transmission systems (FACTS) 
devices are proposed to implement to power grid for 
improving the stability of the systems. In which, a new 
topology structure of a static VAR compensator (SVC), i.e., a 
low-cost improved thyristor-controlled reactors by removing 
series circuit elements, proposed in [1] was used to correct the 
unbalanced problems and improve the reliability of power 
systems. As well as, a hybrid SVC presented in [2] was 
utilized to improve reactive power control capability of wind 
farms fed to an AC grid in case of low-voltage conditions. A 
SVC and a thyristor-controlled series compensator (TCSC) 
were employed in a single-machine infinite-bus system to 
improve transient voltage stability of an asynchronous wind 
farm [3]. A SVC was located at the middle point of the tie line 
to damp out inter-area oscillations occurred in a two-area 
four-machine system [4, 5]. Comparing to the performance of 
a SVC, a static synchronous compensator (STATCOM) gave 
a larger contribution to transient margin and low-voltage 
ride-through capability of squirrel-cage induction 
generator-based wind farms using both calculations and 
simulations [6]. In [7], a STATCOM was connected at the 
point of common coupling to maintain stable voltage and to 
improve the power quality by protecting doubly fed induction 
generator (DFIG)-based wind farm connected to a weak grid 
from going offline during and after disturbances. As proposed 
in [8], a new control strategy using a full-power wind 
permanent-magnet generator (PMG) with a STATCOM to 
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provide reactive power compensation can be used to improve 
stability of transient voltage. In [9], performance of a 
PMG-based wind energy system employing a dynamic 
voltage regulator was compared to one of the system 
employing a STATCOM. It is recommended to use a 
STATCOM in systems with large loads where reactive power 
consumption from the grid could cause serious effects on 
connected loads. 

This paper presents a STATCOM/SVC and load is connected 
to a point of common coupling (PCC). They compensate 
reactive power at the connected bus and maintaining the 
voltage at PCC. The main case study in this paper is to design 
PID controllers for STATCOM and SVC to improve 
oscillations of the studied system. 


n. System Configuration and mathematical models 
Fig. 1 shows the configuration of the studied system. A 
STATCOM/SVC with capacity of 20MVA is connected to 
Bus 5 of a three-machine power system. 



Fig. 1. Configuration of the studied system 


The mathematical models of the synchronous generator (SG), 
STATCOM, SVC, network and there pu differential 
equations are presented in subsection below. 

A. Synchronous Generator Model [10] 

The well-known three-machine nine-bus power system which 
is widely used in power system stability studies was shown in 
Fig. 1. The complete parameters of this system and reference 
frame can be found in [10]. In this case study, the SGi model is 
represented by a two-axis model as block diagram that is 
shown in Fig. 2. In this model, the transient effects are 
accounted for while the transient effects are neglected. 
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The pu differential equations for the i-th synchronous 
generator are described as below. 



loads 


This property is used to obtain the network reduction as 
shown below. With the relation between nodal currents and 
voltages is: 

l sys =Y sys \ (5) 

or 



where 


P(E dl ) = J [-E,n-{x ql -x q )l q ] 

(1) 

■•-[4 

IL 2 I L 3] I s the injected current vector at 

qoi 
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generator nodes, 

p(E -) = —— [- E. + E FDi + (x di - X di )l di ] 


T 

T doi 


v 31 =[v, 

V 2 v 3 ] is the voltage vector at the generator 

p{a ) i ) = — \r rni - [l di E di + I qi E qi - (4. - 4 )l di I q 

]- Di co] (3) 

nodes, 


p(S i )=co i -1 
where: 

(4) 

v 32 = [4 

nodes. 

V 2 , V y J is the voltage vector at the eliminated 


1 (rad) and ‘ (rad/s) represent the rotor angular position 
and angular velocity, respectively; 


the synchronous generator, respectively; 

£ 

FDi (pu) is the field voltage; 

I di (pu) and I q! (pu) are the d- and q-axis currents; 

7dm (s) and qo1 (s) are the d- and q-open-circuit transient 
time constants; 

X X 

di (pu) and qi (pu) are the d- and q- transient reactance; 

X x 

di (pu) and qi (pu) are the d- and q- synchronous 
reactance; 

T 


“ (pu) is the mechanical nominal torque; 
‘ is the damping factor; 


B. Multi-Machine System [10] 

Fig. 3 presents the studied multi-machine system. This 
network has three synchronous generators and three loads. 
Assume that the three loads are represented by three constant 
impedances while three generators are represented by three 
active sources. Therefore, all the nodes have zero injection 
currents except for the generator nodes. 


Ya'vj is the original system admittance matrix. In this case, 
Y sys is the desired reduced matrix and it has the dimensions of 
3x3. 


Expanding the above equations, it yields: 


I 3 = Y 31 V 31 +Y 32 V 32 (7) 

0 = Y 33 V 3I+ Y 34 Y 32 ^ (8) 

From these (7) and (8) equations, V 32 can be eliminated to get: 


I> = (Yj,-Y„YJY„)V„ (9) 

Y = (Y„-Y !! Y i Y !! ) (1Q) 

Because of each SG is modeled under its own qd-atds 
reference frame rotating with its rotor, we can transform from 
generator’s qd -axis rotor reference frame to the common 
QD -axis reference frame as follows: 
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C. STATCOM Model 


( 11 ) 

( 12 ) 


The single line diagram of the employed STATCOM is 
presented in Fig. 4. The relationship between output voltage 
of the STATCOM and its parameters can be written by [6,7]: 
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v dsm = Vdcsta ■k m snr Sln (@PCC + a sm ) (13) 

V q St a = Vdcsta £ m sta ■ Si n (^ PCC + ) (14) 

where 

v and v dsla are the pu q- and d-axis voltages at the output 
terminals of the STATCOM, respectively; 
km sta and a , a are, respectively, the modulation index and 
phase angle of the STATCOM; 

Opcc is the voltage phase angle of the common AC bus; 
V dcsta is the pu DC voltage of the DC capacitor Cm. 


Vpcc 



km sta , a s , a 


Fig. 4. Single-line diagram of the employed STATCOM 

The pu DC voltage, current equations of the DC capacitor C m 
can be described by 


adjusted by injecting the reactive power into system from 
SVC. 



Fig. 7 shows the control block diagram of the studied SVC. 
When the voltage at PCC is lower than the reference voltage, 
the value of B svc of the SVC is positive to inject reactive 
power to the system. When the voltage at PCC is higher than 
the reference voltage, the B svc of the SVC is negative to 
absorb reactive power from the power system. Assuming a 
balanced and fundamental-frequency operation, the 
equivalent B svc of the SVC is a function of the firing angle a 
as shown below [6]. 


■L„ 


(15) 

2a - sin (2a) -nl 2--?M 
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I dcsta = iqsta-k ,n S M cos(0 pcc + a mi )+ i dsta .km sla s\n(d PCC + a su ) (16) 
where 

1 d a ta ' s the pu DC current flowing into the positive terminal 
° f Vdata ; 

R m is the pu equivalent resistance considering the equivalent 
electrical losses of the STATCOM 

i and i dsla are the pu q- and J-axis currents flowing into 
the terminals of the STATCOM, respectively. 

Fig. 5 presents the fundamental control block diagram of the 
STATCOM with the PID controllers. The pu DC voltage 
Vdcsta ' s controlled by a sta while the AC voltage v sla is 


where X L and X c are reactance of reactor and capacitor of 
SVC, respectively. 

Vpcc 

BsVC.max 
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Fig. 7. Control block diagram of the SVC 


varied by changing km sta . 



Fig. 5. Control block diagram of the STATCOM 

D. Static VAR Compenstor (SVC) Model [4] 

In this paper, the proposed SVC is referenced from [4], Its 
single line diagram of the SVC with TCR-FC components 
type was presented in Fig. 6. The voltage at the PCC is 


III. Design PID Controller Of Statcom And Svc Using 
Modal Control Theory 

This subsection describes the design procedure and design 
results of the PID damping controllers for the proposed 
STATCOM and SVC to achieve stability improvement of the 
studied system using pole placement technique. The nonlinear 
system equations are linearized around a selected nominal 
operating point to acquire a set of linearized system equations 
in matrix form of [10] 

/?(x) = AX + BU + VW (18) 

Y = CX + DU (19) 

where X is the state vector, Y is the output vector, U is the 
external or compensated input vector, W is the disturbance 
input vector while A, B, C, and D are all constant matrices of 
appropriate dimensions. 

The eigenvalues of the studied three-machine nine-bus system 
are listed in the second column of Table I. In which, the 
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eigenvalues from modes 1 to 10 relate to the modes of the 
three-machine nine-bus. It can be found by examining the 
system eigenvalues listed in the second column of Table 1 that 
all modes of the system are almost fixed on the complex plane 
and the modes 1 and 2 relating to the rotor angle deviation 
between SG#1 to SG#2 (A<5j 2 ) and the rotor angle deviation 
between SG#1 to SG#3 (A<Sj 3 ), respectively, are changed and 
these modes can be improved by damping controllers. 


Table I. Eigenvalues of the studied system 


Modes 

System 

System with 
STATCOM and 
PID 

System with 
SVC and 
PID 

1 

-0.57 ± 715.69 

-0.8 ± 713.0 

-0.8 +713.0 

2 

-0.72 ± j8.62 

-1.0 ± 78.0 

-1.0 + 78.0 

3 

-5.26 + /7.88 

-5.27 ± /7.89 

-5.26 +78.02 

4 

-5.24 ±./7.94 

-5.28 ± /7.93 

-5.34 +77.96 

5 

-5.37 ±./7.96 

-5.38 + /7.97 

-5.24 + /7.94 

6 

-3.53 +./1.11 

-7.43 + /0.78 

-0.44 + /0.77 

7 

-0.45 +/0.82 

-0.32 + /0.88 

-0.37 +/0.57 

8 

-0.37 + /0.53 

-0.37 +70.65 

-0.27 + /0.20 

9 

-0.24 +/0.03 

-0.18 + /0.15 

-2.03, -4.69 

10 

-0.14, -3.22 

-0.24, -3.23 

-0.41, -2.42 

11 

-6.91 

-1.69, -0.23 

-6.71, -0.14 

12 


-100.02, -99.40 

-26.71, 

-13.71 

13 


-0.33 +70.18 

-7.36 + 
714.25 

14 


-3.71 

-3.23 


The control block diagrams of the proposed STATCOM and 
SVC including the designed PID damping controllers were 
shown in Fig. 5 and Fig. 7, respectively. The PID damping 
controller senses the rotor speed deviation of the SG#1 and 
SG#2 (Ary1 2 ) to generate a damping signal ( V c ) in order to 
improve the damping ratios of modes 1 and 2 of studied 
system. Hence, the output signal in (19) is Y = Awn and U = 
V c is the input vector. The transfer function H(p) of the PID 
controller for the STATCOM/SVC is given by 


H(p) = 


U {p)_ yM 

Y{p) A a n (p) 


1(19) 


where p is the Laplace operator; T is the time constant of 


the wash-out term; K p , K f , and K D are the proportional 
gain, integral gain, and derivative gain of the damping 
controller, respectively. 

For improving the damping of the selected modes, in this 
studied system, two new eigenvalues are selected as (-0.8 ± 
j\ 3.0) and (-1.0 ± j'8.0) for mode 1 and mode 2, respectively. 
The design results of the PID damping controller for the 
STATCOM are given as K P = 32.5, K, = 83.4, K n = 23.4, and 
T w = 0.4 s, and the design results of the PID damping 
controller for the SVC are given as K P = 17.7, K, = 127.5, K n 
= -4.3, and T w = 0.5 s. The new eigenvalues of the studied 
system and the proposed STATCOM and SVC joined with 
the designed PID damping controllers are listed in the third 
and the fourth columns of Table I. In which, the modes of the 
proposed STATCOM/SVC and the designed PID damping 
controllers are from 11-14. It can be clearly observed that 


modes 1 and 2 have been exactly positioned on the desired 
locations on the complex plane. 

IV. Time Domain Simulations Of Nonlinear Model 
In this section, time domain simulations are performed by 
using the nonlinear system model developed in the previous 
section to compare the damping characteristics contributed by 
the proposed STATCOM and SVC joined with the designed 
PID damping controllers on stability improvement of the 
studied system under a three-phase short-circuit fault 
happened at the load C belong to Bus 8 of Fig. 1. 



(a) Rotor angle deviation between SG#1 to SG#2 



(b) Rotor angle deviation between SG#1 to SG#3 


(c) Rotor speed of SG#1 

I --- 


(d) Rotor speed of SG#2 



(e) Rotor speed of SG#3 


Fig. 8. Transient responses of the studied system under a 
three-phase short-circuit fault at load C 

The three-phase short-circuit fault is suddenly happened at t = 
1 s and is cleared after 5 cycles. The simulation results of the 
proposed system using MATLAB/SIMULINK toolbox are 
presented in Fig. 8. In this figure, the comparative transient 
responses of the studied system without FACTS devices 
(black lines), with the SVC and PID damping controller 
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(blue lines), and with the STATCOM and PID damping 
controller (red dotted lines) are respectively plotted. In 
which, the rotor angle derivation of the SG#1 to SG#2 and 
SG#1 to SG#3 are shown in Fig. 8(a) and Fig. 8(b), 
respectively. The rotor speeds of three SGs are also shown in 
Figs. 8(c)-(e). It is clearly observed from the comparative 
transient simulation results that the proposed STATCOM 
joined with the designed PID damping controller can offer 
better damping characteristics than the SVC joined with the 
designed PID damping controller. This is due to the fact that 
the maximum capacitive power generated by a STATCOM 
decreases only linearly with the bus voltage but it drops off as 
square of the bus voltage for a SVC. 

V. Conclusion 

This paper has presented the comparative stability 
improvement studied of a multi-machine system using a 
STATCOM and a SVC. To supply the adequate reactive 
power to the system, PID damping controllers for the 
proposed STATCOM and SVC have been designed by using 
pole placement technique to assign two dominant modes of 
the studied system to the desired locations on the complex 
plane. Comparative time-domain simulations of the studied 
system subject to a three-phase short-circuit fault at the other 
load bus have been systematically performed to demonstrate 
the effectiveness of the proposed STATCOM and SVC joined 
with the designed damping controllers. Since a STATCOM 
provides faster control than an SVC and improved control 
range, it can be concluded from the simulation results that the 
proposed STATCOM joined with the designed controller has 
better damping characteristics to improve the performance of 
the studied multi-machine system under a severe operating 
condition 
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